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ABSTRACT 

Upon formation, degenerate He core white dwarfs are surrounded by a radiative H-rich layer primarily sup- 
ported by ideal gas pressure. In this Letter, we examine the effect of this H-rich layer on mass transfer in 
HeH-C/O double white dwarf binaries that will eventually merge and possibly yield a Type la supernova (SN la) 
in the double detonation scenario. Because its thermal profile and equation of state differ from the underlying 
He core, the H-rich layer is transferred stably onto the C/O white dwarf prior to the He core's tidal disrup- 
tion. We find that this material is ejected from the binary system and sweeps up the surrounding interstellar 
medium hundreds to thousands of years before the SN la. The close match between the resulting circumstel- 
lar medium profiles and values inferred from recent observations of circumstellar absorption in SNe la gives 
further credence to the resurgent double detonation scenario. 

Subject headings: binaries: close — novae, cataclysmic variables — nuclear reactions, nucleosynthesis, 
abundances — supernovae: general — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) are responsible for nearl y half 
of the heavy element production in the Universe (Timmes 
et al. 1995|) and serve as p robes of cosmic acceleration ( Riess, 
et al. |1998t|Perknutter et al. 1999). However, the nature of 



their progenitor systems is still a mystery. While researchers 
agree that SNe la involve the explosions of C/O white dwarfs 
(WDs), three main possibilities for the companion to the ex- 
ploding WD remain in contention: a H-rich donor (ty pically 
referred to as a "single degenerate" system; Whelan & Iben| 



l973j|Nomo to'1982!), another C/O WD ("double degenerate"; 
Iben & Tut ukov 19841 |Webbink|1984|), and a H e WD or He- 



burning star ("double detonation"; Livne|1990 l. Despite con- 
certed effort, no consensus has yet been reached. 

Recent theoretical work (Schwab et al. 20121 IShen et al. 



2012]), in an update of initial studies (e.g., Nomoto & Iben 



T985), suggests that the long term evolution of a double de- 
generate system does not yield a SN la but instead leads to the 
formation of a C -rich giant and then p ossibly to a collapse to 
a neutron star (Saio & Nomoto| |1985|l. Furthermore, a wide 



variety of work, both theoretical (e.g.,|S hen & Bildsten|2007 
Ruiter et al. 2009, Kasen 2010) and observational (e.g., |Di 



5tefano,2010,,Gilfanov & Bogdan,2010 ), suggests that single 
degenerate systems cannot be the dominant progenitor chan- 
nel. While the double detonation scenario has been compara- 



tively less well- studied, recent work ( Fink et al.|2007 2010[ 
Sim et al.|2010| l has resurrected interest in them and spurred 



spurred 

ongomg research. 

One point seemingly in favor of single degenerate progen- 
itors has come from recent observations of circumstellar ma- 
terial (CSM) surrounding 10-30% of SNe la at distances of 
0.1-1 pc from the explosion and velocities of 50- 150 km 
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concluded that the presence of CSM this close to the explo- 
sion implies some form of a H-rich single degenerate progen- 
itor, since several variants of this scenario predict mass loss 
histories that might mimic the observed CSM distributions. 

However, we show in this Letter that H-rich material is also 
ejected hundreds to thousands of years prior to the merger of 
a He WD and a C/O WD, which may lead to a SN la in the 
double detonation scenario. For ambient interstellar medium 
(ISM) densities appropriate for spiral galaxies, this ejecta and 
the swept-up ISM yield a CSM profile that matches observed 
distances, velocities, and column densities. We also find that 
the lower ISM density in elliptical galaxies inhibits the for- 
mation of neutral Na in the CSM, which could explain the 
non-detection of neutral Na absorption in SNe la in ellipticals 
without dust lanes. 

We begin in Section |2] by describing the binary evolution 
that produces a He WD with a thin surface H layer We derive 
the properties of the initially stable H-rich mass transfer from 
the He WD to the C/O WD and perform a time-dependent cal- 
culation of this accretion onto the C/O WD in Section|3] Prior 
to the merger of the He WD with the C/O WD, the H-rich ac- 
cretion yields multiple ejection events akin to classical novae, 
which expel material into the surrounding ISM. We model the 
resulting CSM in Section]?] and conclude in Section[5] 

2. EVOLUTIONARY SCENARIO AND PROGENITOR 
CHARACTERISTICS 

Double detonations may arise in systems with a C/O WD 
primary via stable He mass transfer fro m a He-b urning star 
( |Livne|1990"] l or low mass He WD donor ( |Fink et a l. 2007 ), or 
by unstable mass transfer from a hig her mass He W D ( Guill(>] 
Ichon et al.|20 "T0^ or C/O WD donor ( [Pakmor et al.T2013) . The 
subclass of double detonation progenitors we consider in this 
work are binaries with a relatively high mass 0.3-0.45 Mq 
He WD and a 0.9-1.2 M© C/O WD. This range of He WD 
masses is predicted to yield a He detonation either via an ac- 
cretion stream instability in the l ead-up to the mer ger or dur- 
ing the disru ption of the He WD (Guillochon et al. 2010; Dan] 
et al.||20l"2 l, or possibly during the viscous evolution of the 
merger remnant ( ,Schwab et al.,,2012) . Detonations of either 
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Figure 1. Mass fractions vs. mass below surface of a 0.40 Mq He WD, 10* 
yr (top panel) and 3 X lO' yr (bottom panel) after truncating RGB evolution. 
Solid lines correspond to mass fractions of 'H, 'He, "^He, and '''N as labeled. 
Also shown as a dashed line in both panels is the ratio of ideal gas pressure, 
nkT , to total pressure, Ptotal, where n is the combined nuinber density of ions 
and free electrons, and T is the temperature. 

lower or higher mass C/O WDs underproduce or o verproduce 
^^Ni, respectively, as compared to typical SNe la (Sim et al. 



2010 Ruiter et al. 2013 1. Throughout this Letter, we focus on 
aTiHucial binary with a 0.40 Mq He WD and a 1.0 M© C/O 
WD. Binaries with different components will yield quantita- 
tively, but likely not qualitatively, different results, which we 
defer to future work. 

We construct our fiducial .40 Mc^ He WD with the stel- 
lar evolution code MESAp](Pa xton et al.j201l 'i by truncating 
the evolution of a 1.0 Mq star with an mitial metallicity of 
Z = 0.01 as it ascends the red giant branch (RGB). Chemi- 
cal diffusion is active throughout the calculation. When the 
H-deficient core mass reaches 0.399 Mq, mass loss is turned 
on to rapidly remove the convective H envelope, yielding a 
0.40 Mq He core WD. 

Upon formation. He WDs possess a thin, 2 x 10-'*-2 x 
lO""* Mq H-rich remnant surface layer. This layer evolves 
under the action of cooling, chemical diffusion, and nuclear 
burning, yielding a radiative layer containing essentially pure 
H on top of a degenerate He core polluted with '"^N. We refer 



readers to previous studies (e.g., Al thaus et al.|2001 Kaplan 

tal. 2012 ) for a comprehensive description or He WD evo- 
tion and instead describe our two fiducial He WD models, 
which evolve for 10^ and 3 x 10^ yr after the truncation of the 
RGB prior to the onset of mass transfer to the C/O WD. These 
two timescales are chosen to parametrize our uncertainties in 
the double WD binary's separation following common enve- 
lope evolution. The elapsed time between the formation of 
the He WD and the onset of mass transfer depends strongly 
on this separation and can range from a Hubble time for an 
initial separation of a = 1.3 x 10" cm to 10^ yr for a separa- 
tion of 3.8 X 10'" cm. 

The initial separation and resulting merger timescale are 
critical for determining the mass transfer history once accre- 
tion begins. The H layer on an older He WD will be colder 

^ http://mesa.sourceforge.net/ (version 4589; default settings are used un- 
less otherwise noted) 
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Figure 2. Radius (top panel) and absolute value of the differential radial 
response (bottom panel) vs. transferred mass for our two models. Solid 
(dashed) lines in the bottom panel show positive (negative) values of Dot- 
ted lines in b oth panels show zero-temperature analytic relations from |Nauen-| 
|toi|p72) . 

and thinner, and such a system will have a shorter delay be- 
tween the onset of H-rich mass transfer and the disruption 
of the He WD. The cooling and contraction of the WD, as 
well as a small amount of residual H-burning, also affect the 
abundance structure within the surface layers. Figure[T]shows 
abundance profiles for our two He WDs. Since ^He and '"'N 
act as catalysts for H -burning in classical and recu rrent no- 
vae ( |Townsley (feBUd sten 2004; Shen & Bildsten] [2009 1 



proper calculation of the mass transfer history onto the C/O 
WD, which we describe in the next section, should take into 
account both the changing thermal properties and abundance 
profiles of the aging He WD. 

3. TIME DEPENDENT MASS TRANSFER RATE AND MASS 
EJECTION HISTORY 

When mass transfer begins, the He WD's radiative H layer 
is transferred first. We follow the He WD donor's response 
to this mass transfer by removing mass at a con stant rate of 
10"^ Mq yr"' in MESA (see Kaplan et al.|201 2|for a similar 
analysis of lower mass He WDs^ I'he choice of mass re- 
moval rate in this calculation is essentially arbitrary because 
the thermal timescale in the H layer is > 10^ yr; thus, for 
any accretion rate > 10~'° Mq yr"', the thermal conditions in 
the bulk of the envelope are fixed at the start of mass trans- 
fer. Calculations with constant mass removal rates of 10"^ 
and 10"^ Mq yr"' were also run as a check, with negligible 
differences in the results. 

The donor radius vs. removed mass is shown in the top 
panel of Figure [2] for our two models as labeled. The younger 
model remains hotter and thus has a larger radius. Both mod- 
els have initial radii that are significantly l arger than the zero - 
temperature analytic relation {dotted line; |Nauenberg|| 1972[ ). 
Once the lO""*- 10"~ Mq ideal gas, radiative layer has been 
stripped away, the degenerate He core expands upon further 
mass loss. 

The bottom panel of Figure |2] shows the absolute value of 
the differential radial response, ^ = dXnRjdXnM. It is large 
and positive during the removal of the radiative H layer and 
approaches as this layer is exhausted. As the radius of the 
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Figure 3. Derived mass transfer rate {top panel) and mass ejected from the 
system (bottom panel) vs. time prior to tiie merger for our two models. Each 
system experiences multiple ejection episodes, yielding a total ejected mass 
of 3-6 X 10-5 Mq. 

mass-losing donor expands, ^ becomes negative and eventu- 
ally approaches the value for a zero-temperature low mass 
WD ~ -1 /2), shown as the dotted line. 

In addition to depending on the donor radius and differ- 
ential response, the actual accretion rate must also take into 
account the efficiency of angular momentum exchange be- 
tween the accretor and the orbit (Nelemans et al. 2001 , Marsh 
|et al.|20 04). Assuming conservative mass transfer and angu- 
lar momentum loss only due to gravitational wave radiation, 
the mass transfer rate is 



M2 _ 32G^ MiM2(Mi +M2) 



6 2 



■q-fiq) 



(1) 



The extra term f(q), which accounts for the efficiency o f an- 
gular momentum exchange ([Verbunt & Rappaport|1988| ), has 
little effect on the pre-merger evolution. Eor our binary, the 
derived mass transfer rates for perfectly efficient and ineffi- 
cient angular momentum feedback differ by less than 10% 
while ^ > 10. For ^ = 1, their ratio is 1 .82. Since mass transfer 
in our fiducial binary proceeds at al l times via direct impact 
accretion instead of disk accretion (Lubow & Shu 1975 ), we 
assume that angular momentum feedback is inefficient for the 
remainder of this Letter. 

3.1. Accretion onto the C/0 WD 

Using equation ([TJ, the donor's radial response shown in 
Figure [2] and a relation between the donor's Roche lobe and 
the orbital separation ( |Eggleton|1983 1, we can derive the ac- 
tual mass transfer rate as a function of time. We show this 
accretion rate history in the top panel of Figure |3] for our 
two models. The accretion rate is plotted vs. the time prior 
to the WD merger, which we assume occurs when the accre- 
tion rate reaches the 1.0 Mq C/O WD's Eddington rate of 
2x 10-5 M© yr-10 

^ We assume that super-Eddington accretion will be driven from the system 
via dynamical fricti on with the donor The resulting evolution of the Roche 
overfill factor (e.g., 'Mars h et al.|2004j) implies a tiine of 1 — 10 yr between 
the onset of super-Eddmgton accretion and the donor's disruption. However, 



We model this time-dependent mass transfer onto a 1 .0 Mq 
C/O WD in MESA, taking into account the changing compo- 
sition of accreted material shown in Figure [T[ but ignoring the 
C/O WD's own relatively small H- and He-nch layers for sim- 
plicity. For faster convergence, the optical depth of the outer 
zone is moved inwards to 10^. As H-rich material piles up on 
the surface of the C/O WD, the density and temperature at the 
base of the accreted layer increase until convective H-burning 
is ignited, as in a classical nova. This causes the WD's radius 
to expand until it overflows its Roche radius. 

We assume this material is driven out of the system by its 
interaction with the companio n He WD via a common enve- 
lope (e.g., ILivio et al. 1990 1, and that the ejection velocity 
is roughly equal to the He WD's circular velocity of ~ 1500 
km s"'. The bottom panel of Figure p shows the mass ejec- 
tion history prior to the disruption ofthe He WD. A total of 
3-6 X 10-5 Mq of material is ejected from the binary in mul- 
tiple ejections over the course of 200- 1400 yr prior to the 
merger. The evolution of this ejecta and the swept-up ISM is 
the subject of the next section. 

4. EJECTA - ISM INTERACTION 

The interaction of an expanding shell of material with the 
surrounding ISM has been stu died in detail with respect to SN 
remnants (e.g. Chevalier 1977 ), classical and recurrent novae 
(e.g., Moore & Bildsten 2011), and tidal tail s from double 
white dwarf mergers ( [Raskin & Kasen||2013| l. Analytic so- 
lutions for such evolution are well-known; however, our sit- 
uation is complicated by the existence of multiple ejection 
episodes, which will interact with previously shocked ISM. 
We thus defer discussion of analytic results to future work 
and instead present numerical hydrodynamic simulations. 

We employ a ID, spherically symmetric, Eulerian hydro- 
dynamics code that follows the zone-centered evolution of 
mass density, p, momentum density, pv, and energy density 
pv^ /2 + pu, where u is the specific internal energy. The equa- 
tion of state only allows for ideal gas pressure because the 
medium is optically thin. The mean molecular weight is as- 
sumed to be the solar value, /z = 0.6, for the entire domain, 
which has a 10'^ cm spatial resolution. Fluxes are calculated 
with piecewise upwinded finite differ encing. Artificial viscos- 
ity is i ncluded via the prescription of Tscharnuter & Winkler 
( 1979 1. An optically-thin cooling f unction. A, is inc luded for 
r > 10"* Kby fitting to the results of |Gnat & Sternbe rg (2QQl\. 
The qualitative results in this section have been verified by a 
similar Lagrangian hydrodynamics code. 

We calculate the interaction of our two fiducial binaries' 
ejecta with 10^ K ISM at three mass densities: 0.1, 1, and 



lOnip cm 



where m„ is the proton mass. Ejection episodes 



are initiated by increasing the density in the inner rpert = 
3 X 10'^ cm to a constant value of 'iMi^^/A-Kr^^^, where Mej is 
the ejecta mass, and setting the velocity to Vej = 1500 km/s, 
which is roughly the He WD's circular velocity. The per- 
turbed material's initial radius is chosen such that the ejecta 
is still essentially freely expanding. Each ejection event is 
evolved until the next ejection, as prescribed by the ejection 
history in Figure [3] after which the inner zones are again per- 
turbed in a similar manner. 

In Figure |4j we show the mass density vs. distance from 
the SN at the time of merger for 6 simulations. We find that 

uncertainties such as the efficiency of angular momentum feedback and the 
actual Eddington rate may change this merger timescale. 
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Figure 4. Mass density vs. distance at the time of explosion for our six mod- 
els. Solid lines correspond to the younger (10^ yr) models; dashed lines 
show the older (3 X lO' yr) models. Labels denote the ambient ISM density. 
Thick lines demarcate regions where neutral Na irught exist in the CSM, as 
described in the text. 

when the SN la occurs, the outgoing shock is at a distance 
of 0.1-0.6 pc and the just-shocked ISM has velocities of 
50 - 300 km s"' . These velocities correspond to post-shock 
temperatures of 6 x 10'*-2 x 10^ K. Thus, without substantial 
cooling, there will be a negligible amount of blueshifted, neu- 
tral Na. Significant cooling occurs in only one of our six sim- 
ulations, which has the highest ISM density and the longest 
delay between the first ejection and merger. 

However, there will likely be significant clumping in the 
ejecta (e.g., due to Rayleigh-Taylor instabilities), which is 
unresolvable in our ID simulation. Properly capturing the 
clumping will require multi-dimensional simulations, which 
we defer to future work. Since the cooling timescale, fcooi ^ 
kT /nA, is inversely proportional to density, clumping will en- 
hance cooling in the shocked material and may allow for in- 
creased formation of neutral Na. For now, we approximate 
regions where neutral Na may exist in the CSM by the condi- 
tions r < 10"^ K or fcooi < the simulation age at merger; these 
regions are shown in Figure |4] as thick lines. Note that while 
clumping will alter the CSM s density, it will not significantly 
affect its position or velocity. 

In Figure [5] we show the normalized differential column 
density per velocity bin for regions that might contain neu- 
tral Na, using the same approximation as in Figure |4] The 
three models that contain any blueshifted, possibly neutral, 
Na are shown in the three panels as labeled. This metric ap- 
proximates the neutral Na absorption line profile and shows 
systemic velocity offsets of 50-120 km s~' . For comparison, 
the red dashed line in the middle panel shows the normalized 
average of the Na Di and D2 a bsorption profiles f or SN 20071e 
84 d after B-band maximum ( Simon et al.|2009] l. 

The velocity profiles of the possibly neutral Na resem- 
ble those seen in observations, and the Na column densities 
(4 X 10" - 4 X 10'^ cm"-) overlap with those derived from 
SNe la showing variable absorption lines. However, our col- 
umn de nsities are a factor of a few higher than the observed 

Given our ad hoc inclusion of 




mean of Sternberg et al. (201 1 1 
the effects of clumping and the unconsidered complications of 
time-dependent photoionization and recombination following 
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Figure 5. Differential column density per unit velocity, da/dv, vs. velocity 
of Na that might be neutral, as approximated in Figurep] The normalization is 
such that the largest trough with non-zero velocity has a value of 1 . Thermal 
Doppler broadening is not included. The velocity axis is reversed to match 
observational convention. Values for the age of the He WD and the ambient 
ISM density are as labeled. The red dashed line in the middle panel shows an 
average of SN 20071e's Na D absorption profiles 84 d after B-band maximum. 

the SN, our derived neutral Na column densities are merely 
suggestive. Future multi-dimensional work will enable more 
accurate predictions and possibly correct this mismatch. 

No narrow Na absorption has yet been detected in a SN 
la in an elU ptical galaxy without obvious dust lanes ( |Foley| 
et al.||20l"2| . Our results agree with this finding for several 
reasons. Elliptical galaxies have lower ISM densities, which 
yield a smaller column of shocked material at the time of the 
explosion. Furthermore, the lower ISM density implies both 
lower post-shock densities as well as less shock deceleration, 
which means higher CSM temperatures when the SN la oc- 
curs. These factors increase the cooling timescale and de- 
crease the amount of neutral Na. 

5. CONCLUSIONS 

In this Letter, we have considered the effect of the H-rich 
layer that surrounds a He WD on the He WD's interaction 
with a C/O WD companion prior to a SN la. We have calcu- 
lated its structure (Section 2), its impact on mass transfer and 
its ejection (Section 3), and the ejecta's interaction with the 
surrounding ISM (Section 4). We have found that if a SN la 
occurs when the He and C/O WDs merge, the characteristics 
of the CSM at the time of the explosion match recent observa- 
tions of neutral Na surrounding 10-30% of SNe la in spiral 
galaxies. We have also found that the lower ISM density in 
elliptical galaxies inhibits the formation of significant neutral 
Na in the CSM, which may be why these features have not 
been detected in such SNe la. 

SN ejecta have been obse rved to colhde with surrounding 
CSM in several SNe la (e.g. , |Hamuy et al.|2003||Dilday et al. 
|2012[ [Silverman et al.|2013] ). ITiis interaction requires signif- 
icant material within < 10^'' cm, which is difficult to produce 
in our model unless the SN la occurs < 2 yr after an ejection 
event. More detailed modeling and inclusion of the super- 
Eddington accretion phase just before the merger may help to 
shed light on these observations. 

While our results are promising, future studies are neces- 
sary to strengthen the findings. Further work will include an 
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exploration of a range of WD masses and ages, simulations 
of the ejecta - ISM interaction in multiple dimensions, which 
will allow for clumping and non-spherical ejection, and cal- 
culations of the time-dependent photoionization and recombi- 
nation after the SN la's UV flash. The effects of tidal heating, 
which will likely be si gnificant for these extremely close bi- 
naries (e.g., Piro|201 1 1, should also be considered. 
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